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' ^'Abstract 

the calorimetric determination of the primary energy of extensive air showers, measured by fluorescence telescopes, a precise 
^ knowledge of the conversion factor (fluorescence yield) between the deposited energy in the atmosphere and the number of emitted 
5^ fluorescence photons is essential. The fluorescence yield depends on the pressure and the temperature of the air as well as on the 

Water vapor concentration. Within the scope of this work the fluorescence yield for the eight strongest nitrogen emission bands 
Ci between 300 nm and 400 nm has been measured using electrons from a ^°Sr-source with energies between 250 keV and 2000 keV. 
' 'Measurements have been performed in dry air, pure nitrogen, and a nitrogen-oxygen mixture at pressures ranging from 2 hPa to 

990 hPa. Furthermore the influence of water vapor has been studied. A new approach for the parametrization of the fluorescence 
^ yield was used to analyze the data, leading to a consistent description of the fluorescence yield with a minimal set of parameters. 

The resulting absolute accuracies for the single nitrogen bands are in the order of 15 %. In the investigated energy range, the 

fluorescence yield proved to be independent of the energy of the ionizing electrons. 
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1. Introduction 

, The measurement of air fluorescence is used by many 
modern experiments (i.e. Pierre Auger Observatory [1 , 
lliRes [2]) to detect extensive air showers (EAS), induced 
by ultra-high energy cosmic rays. The secondary EAS par- 
ticles, predominantly electrons and positrons, deposit their 
energy in the atmosphere by exciting or ionizing the air 
molecules which afterwards may relax by emitting fluores- 
cence photons. As pointed out by Bunner 3 most of these 
emissions, in the wavelength range between 300 nm and 
400 nm, originate from transitions of neutral or ionized ni- 
trogen molecules. These faint emissions can be measured 
by fluorescence telescopes, allowing the observation of the 
longitudinal development of EAS through the atmosphere 
and a calorimetric determination of the primary cosmic- 
ray energy. The conversion factor between the deposited 
energy in the air and the number of emitted fluorescence 
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photons is the so-called fluorescence yield Y\{p,T) which 
depends on the air pressure p and temperature T at the 
place of emission, as well as on the wavelength A of the emit- 
ted photons. If Y\{p, T) does not depend on the energy of 
the exciting particles the number of fluorescence photons, 
observed by a telescope, is directly related to the deposited 
energy according to 



dx 



dEde 

dx 



A 



Yx{p,T) -Tatmi^x) ■ edet{\x) (1) 



where dx denotes an interval along the shower axis, 
TatmCAjCc) is the atmospheric transmission factor and 
£det(A,a;) corresponds to the detection efficiency of the 
telescope. The summation goes over all wavelengths A 
passing the filter of the telescope, typically ranging from 
300 nm to 400 nm. The one dimensional representation 
above only holds for large distances between the EAS and 
the telescope where the lateral spread of the EAS can be 
neglected. Integrating Eq. ([T]) over x results in the total 
electromagnetic energy of the EAS and thus is the most 
direct measure of the primary cosmic ray energy. The pre- 
cision of this method is, however, limited by the present 
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Fig. 1. Scheme of the energy states of the 2P and IN electronic-vi- 
brational band systems of N2 and N J . The corresponding spectral 
bands in Fig. [2] are drawn in the same colors. 

uncertainties of the fluorescence yield of about 15 % to 
30 % and the lack of knowledge about its energy depen- 
dence. In recent years this gave rise to a number of new 
laboratory experiments, i.e. Kakimoto et al. [4], Nagano et 
al. AIRFLY T, FLASH fg; or AirLight, aiming at a 
precise measurement of the fluorescence yield over a wide 
energy, pressure and temperature range. 

This paper reports on the data analysis and the results 
of the AirLight experiment at Forschungszentrum Karl- 
sruhe in Germany. Section [5] describes the process of nitro- 
gen fluorescence in air and introduces a new approach to 
model the fluorescence yield in a consistent way and with a 
minimal set of parameters. The experimental setup is ad- 
dressed in Section [3] and finally the data analysis and the 
results are discussed in Section |4] and [S] This paper is ex- 
tracted from the Ph.D. thesis [S] of the corresponding au- 
thor which can be referred for more details. Some changes 
have been applied to correct minor errors and to improve 
the readability. 

2. Nitrogen Fluorescence in Air 

Nearly all the air fluorescence emissions in the wave- 
length range between 300 nm and 400 nm originate from 
neutral or ionized nitrogen molecules [3111] . The fluores- 
cence spectrum of molecular nitrogen is a band spectrum. 
In contrast to atomic line spectra, molecular spectra consist 
of a variety of broad bands. This band structure is caused by 
the vibrational and rotational movements of the molecular 
nuclei which modify the energy states of the electrons. For 
diatomic, homo-nuclear molecules the energy of a molecu- 
lar state can be expressed as the sum of three contributions 
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where E^i is the potential electron energy of a static 
molecule and Eyib and Erot are additional contributions 
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Fig. 2. Nitrogen fluorescence spectrum between 300 nm and 400 nm 
in dry air at 1013 hPa measured by Ulrich et al. |10) . 

due to vibrations and rotations of the molecular nuclei. 
This approach is known as the Born-Oppenheimer approxi- 
mation |12|13j . The energy scales of the three contributions 
are very different and approximately behave like 



Ep] : E, 



vib 



Erot — 1 




(3) 



where m is the electron mass and M denotes the mass of 
the molecular nuclei [12]. This causes an electronic molec- 
ular state to split into several vibrational levels and every 
vibrational level again has a rotational substructure, as is 
illustrated in Fig.[T] Very often the rotational substructure 
is not resolved, which leads to the observation of vibra- 
tional bands with a sharp edge at one side and a shading 
to the other side, as can be seen in Fig. [2l Since the ro- 
tational substructure plays no role for the application on 
EAS measurements it will be disregarded in the following. 

An "electronic band system" is established by all 
electronic-vibrational transitions having the same initial 
and final electronic states. The nitrogen fluorescence spec- 
trum between 300 nm and 400 nm nearly entirely consists 
of transitions of the "second positive" (2P) band sys- 
tcn{3 of N2 and the "first negative" (IN) band system of 
^2 i3>llj. The spectroscopic notation for these transitions 
is given by 

- 2P-System: C^Uuiv') B^Ug{v") 

- IN-Systcm: B^Y.+ {v') -> X^j:j{v") 

where the quantum numbers v' and v" of the initial and 
final vibrational levels are written within the parentheses 
behind the symbols for the electron configuration of the 
corresponding molecular state. Accordingly these transi- 



^ The notation "first negative" and "second positive" system indi- 
cates the place of appearance of the corresponding light emissions in 
gas discharge tubes. "Negative" systems appear closer to the cath- 
ode whereas "positive" systems are attracted by the anode. 



tions are denoted as 2P{v',v") in case of the 2P-System 
and lN('i;', v") for the IN-System. In the following the sym- 
bols u' and v" always denote initial and final electronic- 
vibrational states of a certain electronic band system. Any 
numerical values of v' or v" indicate the vibrational level of 
the corresponding electronic state. Pure vibrational tran- 
sitions are not a subject of the following discussions. 

2.1. De-Excitation 

Assuming a number of excited nitrogen molecules in a 
certain electronic-vibrational state v' and neglecting any 
kind of cascading into this state, the deactivation rate can 
be calculated according to the general decay law 

^ = -(A„.t + 5]A,',.") -^.'W , (4) 



where Nyi (t) denotes the number of molecules in the state 
w' at a certain time t. The total decay constant A^', which 
equals the reciprocal mean lifetime Ty> of the state v' , is 
the sum of the decay constants of all potential decay chan- 
nels. The quantities Ayiy denote the decay constants or 
transition probabilities for all radiative transitions v' v" 
within the underlying band system and are also known as 
Einstein coefficients. A detailed list of Einstein coefficients 
of many nitrogen and oxygen band systems can be found in 
the comprehensive work of Gilmore et al. [14j . In addition, 
there might be any kind of other radiating or non-radiating 
transitions which are summarized in the common decay 
constant Ar-est- From Eq. ([4]) two fundamental properties 
can be immediately derived: 

(i) The measurable mean lifetime r„' of an electronic- 
vibrational state v' is the same for all transitions 
v' —f v" originating from this state. 

(ii) The intensity ratios between spontaneous transitions 
of the same electronic-vibrational state v' are con- 
stant and correspond to the ratios of their Einstein 
coefficients Ayi^^jn. 

The latter property causes the 2P and IN band systems of 
the nitrogen spectrum to be assembled from several sub- 
spectra for each vibrational level w', as is illustrated in 
Fig. [21 In the wavelength range between 300 nm and 400 nm 
the most intensive nitrogen transitions of the 2P system 
originate from the vibrational levels v' — 0,1, whereas 
the IN system has just one notable transition 1N(0,0) at 
391.4 nm. For each of these sub-spectra the intensity ratios 
between their vibrational bands are constant, but the abso- 
lute normalization of each sub-spectrum in general changes 
differently with pressure and temperature, as will be ex- 
plained in the following. 



2.2. Radiationless Deactivations: Quenching 

An excited electronic-vibrational state v' can also be- 
come de-activated by radiationless processes such as ro- 
tational, vibrational or translational energy transfer dur- 
ing collisions with other molecules fT3\. These so-called 
quenching processes strongly depend on the number den- 
sity (— > Pressure) and the velocity (— > Temperature) of the 
colliding molecules. 

Assume an excited nitrogen molecule in the state v' hit- 
ting a gas target where the target molecules are at rest. 
Under such conditions a coUisional deactivation constant 
A" J, can be defined according to 

K[x^n^-'^N2A^rel)-Vrel , (5) 

where Ux is the number density of the target molecules of 
type X and a'^^ xi'^^rei) denotes the total coUisional cross- 
section for the deactivation of a nitrogen molecule in the 
state v' which was injected with a relative velocity Vrei with 
respect to the target molecules. In reality, the projectile is 
part of the target gas and the Maxwell-Boltzmann distri- 
bution applies for the velocities of the target and projectile 
molecules. Therefore Eq. ([5]) has to be averaged over all rel- 
ative velocities Vrei between the colliding molecules j3l9| . 
The coUisional cross-sections xi''^rei) in general depend 
on the relative velocities Vrei and thus on the temperature 
of the gas. At present, no temperature dependence has been 
reported and therefore ^ in the following is assumed 
to be constant. In this approximation, averaging over the 
molecular velocities in Eq. ^ leads to the expression 

K,x = nx ■ cr^^^x ■ {VN2) ■ \ , (6) 

V nijx 

V ' 

= : QS'(T) 



where (WJV2) — y^SkT/irmpf^ is the mean velocity of the 
nitrogen molecules at temperature T according to the 
Maxwell-Boltzmann distribution and rux is the mass of the 
target molecules of type x. It is convenient to introduce a 
quenching rate constant Qj; (T) for each gas constituent 
x. Since (T) depends on (vn^) it is not a real constant, 
but it is proportional to Vt if no additional temperature 
dependence from the coUisional cross-sections aj^ ^ has to 
be taken into account. All the quenching rate constants in 
this paper are quoted for a temperature of 293 K (~ 20° C). 
Quenching rate constants for other temperatures can be 
computed according to 

QtiT)^/^^-Q^mK) . (7) 

If the gas is a mixture of different constituents x, the to- 
tal coUisional deactivation constant A" is the sum of the 
deactivation constants of all gas constituents: 



Q 



X X 

The last expression is proportional to p/ y/T by a factor Cv' 
and was obtained by substituting the number densities 
through the partial pressures px = n^kT, as follows from 
the ideal gas law. In a second step the partial pressures Px 
have been expressed as Px = fx ■ where denotes the 
relative abundance of the different gas constituents, e.g. for 
dry air /n, = 0.78, fo, = 0.21 and /av = 0-01. 

With this definition for the coUisional deactivation con- 
stant A" {p, T) the total decay constant or reciprocal life- 
time 1 /r^' of an electronic vibrational state v' can be ex- 
pressed as 

— ^ = A^'(p,r)-)-A;;:,, + ^A„,,„. , (9) 



where A"g^j incorporates the decay constants for all re- 
maining deactivations or transitions with constant transi- 
tion probabilities. At zero pressure A^ {p, T) vanishes and 
the total lifetime t^' (p, T)\^^q becomes equal to the intrin- 
sic lifetime Tq , which is the reciprocal sum of all constant 
transition probabilities. 

Instead of the quenching rate constants {T), it is 
sometimes more convenient to introduce a reference pres- 
sure p'^, (T) for a given gas mixture. The reference pressure 
Pyi{T) is defined as the pressure where the collisional de- 
activation constant AJ! (p, T) equals the reciprocal intrinsic 
lifetime I/tq . With this definition Eq. ^ transforms to 
the simple form 

which immediately reveals the linear pressure dependence 
of the reciprocal lifetime 1 /t„/ {p, T) for a constant temper- 
ature. 

2.3. Excitation: Intrinsic Fluorescence Yield 

The excitation of the relevant molecular nitrogen states 
in EAS generally occurs via many different processes. The 
initial electronic state of the IN system is usually 

activated through the direct ionization of neutral nitrogen 
molecules by high energy electrons 

e + N2(xiS]+) — >N2+*(B2s+)-|-e + e , (11) 

whereas the excitation of the C^Iiu state of the 2P system 
mainly occurs via low energy secondary processes such as 
the recombination of ionized nitrogen molecules 



e + N+ (X2S+) N2* {C^Uu) , (12) 

or the direct excitation from the ground state 

e(t) + N2(Xil]+) ^N2*(C3n„)-|-e(i) . (13) 

The latter process, to the first order, is forbidden, since it 
requires a change in the spin quantum number of the molec- 
ular state. Nevertheless, this excitation becomes possible 
for low energy electrons via electron exchange if spin-orbit- 
coupling is taken into account [3I13| . 

In general, an electron releases its energy in a series of 
different processes and only a small fraction of the en- 
ergy, which was deposited in a certain volume, is finally 
converted into fluorescence photons. Thus if H„' {E) is the 
number of nitrogen molecules per deposited energy which 
have been excited into the state u', the fluorescence yield 
Yyi {E,p, T) for the radiative transition v' — > v" , in units 
of photons per deposited energy, can be expressed as 

Y,,y,iE,p,T) = K'{E) ■ A,,y, ■ T,.{p,T) 

= ^,,{E)-A,,y,-r^' -^^^^^fl . (14) 
" « ' '^0 

=■■ K',."<^^^ 

The term Ayi ■ Tq in the flnal expression determines the 
probability for an excited state v' to relax via the transi- 
tion v' v" if collisional quenching effects are neglected. 
Therefore, the intrinsic fluorescence yield Y^, yii{E) is de- 
fined as the number of emitted fiuorescence photons per de- 
posited energy in the absence of collisional quenching. The 
collisional quenching is finally taken into account by the 
separate factor Tyi (p, T) / Tq which determines the fraction 
of excited states v' which are not going to be deactivated 
by collisional processes. This approach naturally separates 
between excitation and de-excitation processes. Potential 
energy dependencies of the effective excitation process only 
affect the intrinsic fluorescence yield Y^, ^„{E), whereas 
pressure and temperature only act on the quenching and 
thus on the lifetime Tyi {p, T) . 

As pointed out in Section [2Tl the intensity ratios between 
transitions emerging from the same electronic-vibrational 
state v' are always constant. In order to apply this relation 
to Eq. (fT4|) , a main-transition v' ^ k needs to be specified 
for each initial vibrational state v' . A main- transition is de- 
fined as the most intensive transition of the corresponding 
vibrational sub-system, but in principle the choice is com- 
pletely arbitrary. In this work the main-transitions for the 
investigated band systems have been chosen to be 2P(0,0), 
2P(1,0) and 1N(0,0). Now the intensity ratio Ry'y of a 
transition v' — )■ v" can be quoted with respect to the 
main-transition as the ratio of the corresponding Einstein- 
coefficients: 

Rv' ,v" — ~1 J (1^) 
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Fig. 3. Sketch of the AirLight chamber. After exiting the coUimator 
the electrons traverse 10 cm of gas before they are stopped in the 
scintillator. Seven PMTs, equipped with different filters, measure the 
fluorescence emissions perpendicularly to the electron beam. 

All main-transitions have Rv',k = 1- With these definitions 
Eq. (|14p transforms to the final expression 

Y,,y,{E,p,T)^Y,\^{E)-R,,y,-^^^^^^ , (16) 

^0 

where Y^, j^{E) denotes the intrinsic fiuorescence yield of 
the main-transition. All quantities in this equation now 
have a clear physical meaning and are able to be measured 
experimentally. Furthermore, all transitions emerging from 
the same initial state v' are treated consistently with a 
minimal set of parameters. 

In previous papers the fluorescence yield is often defined 
as the number of emitted fluorescence photons per meter 
track length of the ionizing particle. This definition has 
several disadvantages as will be explained in the following. 
First of all it requires all secondary electrons to be tracked 
as well in order to obtain correct results. Since the field of 
view of all laboratory experiments is limited to a relatively 
small region around the electron beam, secondary electrons 
can easily escape and thus do not contribute to the fiuo- 
rescence emissions in this volume. It is even possible that 
electrons are scattered back from the chamber walls into 
the field of view and induce additional fluorescence pho- 
tons. All these effects are very difficult to treat if the fiuo- 
rescence yield is defined as photons per meter. Instead it is 
relatively simple to retrieve the deposited energy in a cer- 
tain volume by Monte Carlo simulations (see Section l4.3.ip 
and to determine the number of fluorescence photons per 
deposited energy, which accounts for all effects mentioned 
above. Another drawback of the definition in photons per 
meter is that the fluorescence yield depends on the under- 
lying energy loss model as well as on the energy and the 
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Fig. 4. Effective transmission curves of the interference filters (solid 
lines) and the broad band M-UG6 filter (dashed line) overlaid with 
a nitrogen fluorescence spectrum measured by Ulrich et al. |10j . 

pressure, which makes it rather difficult to compare values 
of different authors. All these problems diminish if the flu- 
orescence yield is defined as photons per deposited energy, 
which seems also to be more natural, since the deposited 
energy is the maximum energy which may dissipate into 
fluorescence photons. This definition is much more conve- 
nient, particularly for EAS measurements, since it allows 
the number of observed photons to be directly converted 
into an energy deposit profile of the EAS by means of Eq.[TJ 
Therefore, the fiuorescence yield in this paper is always de- 
fined as photons per deposited energy. The " old" definition 
in terms of photons per meter can be obtained by multiply- 
ing the fluorescence yield with an appropriate energy loss 
function 

3. The AirLight Experiment 

The AirLight experiment |9j at Forschungszentrum Karl- 
sruhe was designed to measure all the different quantities 
in Eq. (|16p . The experimental setup is similar to the exper- 
iments done by Kakimoto and Nagano ct al. [41516] . As is 
shown in Fig.[3]it consists of a cylindrical aluminum cham- 
ber in which electrons are injected along the chamber axis. 
The electrons are emitted from a ^''Sr-source situated at the 
top of the chamber. The source has an activity of 37 MBq 
with an end point energy of 2.3 MeV. The electrons are 
collimated by massive lead rings with an inner diameter of 
5 mm, before they enter the inner volume of the chamber. 
After traversing 10 cm of gas (dry air, pure nitrogen, or 
a nitrogen-oxygen mixture) they are finally stopped in a 
plastic scintillator of 2.5 cm height and 4 cm diameter [S]. 
The plastic scintillator measures the energy of the electrons 
with an energy resolution of about 10 % at 1 MeV. The 
electron rate at the scintillator varies between 10 kHz and 
20 kHz, depending on the pressure in the chamber which 
can be adjusted between 2 hPa and 1000 hPa. 

Due to the relatively low electron rate the intensity of 
the electron-induced fluorescence light is not sufficient to 
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Fig. 5. Exponential time distribution between the electron signals 
in the scintillator and the photon signals in channel 3. The equally 
distributed offset is due to accidental coincidences (thermal noise) 
in the PMT. The total fluorescence signal in this channel is a super- 
position of the 2P(0,2) and the 2P(1,3) transitions. 
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Fig. 6. Single photoelectron distribution of channel 3. Only signals 
between 0.5 pe and 2 pe are considered in the data analysis. Roughly 
11.6 % of the signals are between the discriminator threshold and 
0.5 pe (yellow region). The gray region is a crude estimate of 15 % 
for the signals below the discriminator threshold. 



allow the use of a spectrometer or monochromator to make 
a spectral resolved measurement. Therefore, seven 2" Pho- 
tonis photomultiplier^ (PMT) equipped with different fil- 
ters are mounted perpendicularly, at a distance of 20 cm, 
to the chamber axis to measure the fluorescence photons in 
several wavelength ranges. The active area of each of the 
photocathodes is defined by a 4 cm aperture. As shown in 
Fig.m the set of filters consists of one broad band M-UG6 
absorption filter ranging from 300 nm to 410 nm, identical 
to those used in the fiuorescence telescopes of the Pierre 
Auger Observatory [T] , and six narrow band interference fil- 
ters matched to the most prominent nitrogen bands. Since 
interference filters are known to change their transmission 
characteristics with the incident angle of the light, the ef- 
fective transmission functions plotted in Fig. [?] were ob- 
tained by averaging over the angle distribution of the fiuo- 
rescence photons [1519] which was determined by means of 
Monte Carlo simulations f9]. The pressure dependence of 
the angular distribution, due to the multiple scattering of 
the electrons, turned out to be negligible. 

The experiment measures coincidences between individ- 
ual electron signals in the scintillator and single photon 
signals in any of the PMTs. The coincidence condition is 
fulfilled if a photon signal appears within an interval of 
120 ns after an electron signal was detected. The distribu- 
tion of time differences (see Fig. [5]) between electron and 
photon signals, as well as the single photoelectron distribu- 
tions (see Fig. [6]), are sampled for each filter channel by a 
CAEN V488 TDC[I] and a Lecroy 1182 ADCH- One ADC 
channel is sampling the energy distributions of the elec- 
trons. After every tenth coincidence, the coincidence con- 
dition is disabled for the next electron to measure also the 
unbiased (free) electron energy spectrum. Finally the ab- 
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solute numbers of free and coincident photon and electron 
triggers are counted by a CAEN V260 scaler module with 
a guaranteed input frequency of 100 MHz. The typical free 
rates of the photon detectors at 20 °C ranged from 400 Hz 
to 600 Hz. 

By means of the TDC time spectra, fluorescence signals 
can be clearly discriminated from accidental coincidences. 
As is shown in Fig. [5] the time distribution of the fiuores- 
cence photons follows an exponential decay law as is ex- 
pected from Eq. (|4]). Contrary to this, accidental coinci- 
dences show a flat distribution, since the timing of the elec- 
tron and photon signals is completely uncorrelated. The 
time resolution at of each filter channel is of the order of 
0.7 ns. Therefore, the time distribution of the fiuorescence 
signals can be described by a convolution of an exponential- 
and a Gaussian function with a = ut which is given by the 
expression 

dN 17V ^ ( t^-t+^ \ 

—— ~ - — • e ^ • • eric ■= — . (17) 

dt 2t y V2(7t j ^ ' 

The time integral of this function corresponds to the total 
number of fluorescence photons N , whereas r denotes the 
lifetime of the emitting nitrogen state. As can be seen in 
Fig.[5]the time origin to is in the order of 16 ns and is roughly 
the same in all filter channels. In general there are several 
nitrogen bands with different lifetimes r overlapping in one 
filter channel. Therefore, the measured time distribution 
results from a superposition of several distributions given 
by Eq. (|17p with different lifetimes and intensities but with 
the same values for the time origin and time resolution of 
the relevant filter channel. 

The number of detected photons N\^^ in a filter channel i 
is linked to the real number of emitted photons N through 
the relation 
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The integral in this expression corresponds to the spec- 
tral efficiency , which contains the nominal quantum effi- 
ciency function £q£;(A) of the PMT, as quoted by the man- 
ufacturer, the effective transmission function T(A) of the 
filter and the normalized spectral distribution ^ of the 
emitted photons. If the spectral distribution is known, the 
integral e* can be easily computed, since the unknown nor- 
malization factor fi of the individual quantum efficiency 
function eQ^{X) = fi ■ eQ^{X) has been moved out of the 
integral. The PMT acceptance is given by the geometri- 
cal factor £f2 which is the same for all filter channels, since 
the whole setup is completely symmetric against the cham- 
ber axis. For fluorescence measurements en has been deter- 
mined by means of Monte Carlo simulations to a value of 
(0.2586 ± 0.0003) % which doesn't depend on the pressure. 
The collection efficiency e^^i of the PMTs is assumed to be 
constant over the relatively small filter range. Finally we 
must account for the cut efficienc y g' ,t, since only photon 
signals between 0.5 pe and 2.0 pq3 are used in the data 
analysis as illustrated in Fig. [51 The normalization factor 
fi, the collection efficiency e*^; and the cut efficiency are 
the only unknown quantities in Eq. (|18p . For convenience 
they have been merged into the final calibration constant 

fi 

J cal ' 



3.1. Absolute Calibration 

In principle the calibration constants /*^; can be indi- 
vidually determined by a calibration measurement where 
a well calibrated light source is placed in the center of the 
chamber. In the present work it was not possible to per- 
form such an end-to-end calibration since an appropriate 
light source was not yet available. Nevertheless a relative 
calibration has been performed where an uncalibrated but 
stable and axial symmetric light source with a well known 
spectrum was used instead. Furthermore the filters in front 
of the PMTs have been removed, resulting in T(A) = 1 for 
each channel. This way all PMTs are illuminated with an 
equal but unknown amount of photons and the acceptance 
en is the same for all channels due to symmetry reasons. 
By means of Eq. ([TS]) the calibration constants fl^i have 
been relatively determined according to 



the fraction of the detected photons with respect to the 
reference channel times the reference constant. In order to 
evaluate absolute values for the calibration constants, the 
reference constant /^^^ was estimated under the following 
conditions: 

The shape of the single photoelectron distribution of the 
reference channel, shown in Fig. [51 as well as the discrimi- 
nator threshold are assumed to be stable over the whole pe- 
riod of measurements. Stable photoelectron distributions 
are also required for all the other channels. Furthermore, 
the photoelectron distributions of accidental coincidences, 
which are mainly caused by thermal noise of the PMTs, 
must be roughly the same as for real photon signals. These 
requirements have been proven to be satisfactorily fulfilled 
over the whole data-taking period [9] . 

The area between 0.5 pe and 2.0 pe in the single photo- 
electron distribution shown in Fig. [51 corresponds to 100 % 
of the detected photons which have been considered in the 
data-analysis. Signals above 2.0 pe are mainly due to noise 
or background and can be neglected. However, roughly 12 % 
of the signals are situated between the discriminator thresh- 
old and the quality cut at 0.5 pe (yellow region) and have 
to be taken into account. The number of signals between 
the pedestal at pe and the discriminator threshold (gray 
region) is unknown but has been conservatively estimated 
to be less than 15 % of the detected signals. Therefore, half 
of this value (7.5 %) is added to the number of detected sig- 
nals and the other 7.5 % are taken as the systematic uncer- 
tainty of this estimation. Thus, if Ncut denotes the number 
of signals ranging from 0.5 pe to 2.0 pe, the real number of 
signals Nreai can be evaluated by the expression 



Nreal = (100 % +12% + 7.5 % ± 7.5 %) • N^ut 



(20) 



which finally results in an estimation for the cut efficiency 
of channel 3: 



1 



Nreal (1 . 195 ± 0.075) 



0.837 ±0.053 



(21) 



The PMT quantum efficiency function Sq^{X) of chan- 
nel 3, with a maximum value of 26.6 % at 420 nm, is as- 
sumed to be known with an uncertainty of 5 %, which im- 
plies a normalization factor f^ = 1.00 ± 0.05. For the col- 
lection efficiency e^^j a value of 0.880 ± 0.088 has been as- 
sumed. Using these values the absolute calibration factor 
fcai reference channel 3 can be estimated according 

to 



n: 



det 



f 



cal 



det 



f3 

•> cal 



J cal 



det 



fcal I (19) 



det 



where channel 3 was chosen for technical reasons to be the 
reference channel. In case of equal PMT types the spectral 
efficiencies e* cancel and the calibration constant is just 



^ [pe]; photoelectrons 



f3 ^ p3 
J cal ^c 



^col 



/a ^ 0.737 ±0.095 



(22) 



with a systematic uncertainty of about 13 %. By means of 
Eq. the calibration constants of the other channels can 
be derived relative to Z^^;. Relative calibration measure- 
ments have been performed once a month over the whole 
data-taking period and turned out to be extremely sta- 
ble [5]. The resulting detection efficiencies e^et = fcai'Sn-Ss 
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Table 1 

Overview of the different filter channels and the corresponding nitrogen bands. The detection efficiencies are determined according to 
^det = fcal ■ ef2 ■ Esi Using a common acceptance value of = 0.2586 %. The last column of the table indicates the relative systematic 
errors Asd^t- 



Channel 


Filter 


fcal 


Band 


A [nm] 


Sa 


Sdet [%o] 


Aedet [%] 





M-UG6 


0.832 ±0.109 


2P(1,0) 


315.9 


0.077 ± 0.006 


0.166 ±0.025 


14.93 








2P(0,0) 


337.1 


0.158 ±0.008 


0.340 ± 0.048 


14.09 








2P(1,2) 


353.7 


0.192 ±0.009 


0.413 ±0.057 


13.86 








2P(0,1) 


357.7 


0.196 ±0.009 


0.422 ± 0.058 


13.84 








2P(1,3) 


375.5 


0.200 ± 0.008 


0.430 ± 0.059 


13.76 








2P(0,2) 


380.5 


0.194 ±0.008 


0.418 ±0.057 


13.74 








1N(0,0) 


391.4 


0.152 ±0.006 


0.327 ±0.045 


13.68 








2P(1,4) 


399.8 


0.085 ± 0.003 


0.183 ±0.025 


13.68 


1 


(a) 


0.714 ±0.094 


2P(1,0) 


315.9 


0.094 ± 0.004 


0.173 ±0.024 


13.99 


2 


(c) 


0.825 ±0.108 


2P(1,2) 


353.7 


0.014 ±0.001 


0.029 ± 0.004 


15.39 








2P(0,1) 


357.7 


0.053 ±0.003 


0.112 ±0.016 


14.03 


3* 


(d) 


0.737 ±0.095 


2P(1,3) 


375.5 


0.078 ± 0.004 


0.148 ±0.020 


13.61 








2P(0,2) 


380.5 


0.106 ±0.005 


0.202 ±0.027 


13.53 


4 


(b) 


0.854 ±0.114 


2P(0,0) 


337.1 


0.060 ± 0.003 


0.133 ±0.019 


14.38 


5 


(e) 


0.788 ±0.104 


1N(0,0) 


391.4 


0.081 ±0.004 


0.165 ±0.023 


13.84 








2P(1,4) 


399.8 


0.042 ± 0.002 


0.085 ±0.012 


14.08 



* Reference channel. 



out before each series of runs with a certain gas. 

4.1. Global Minimization Function 

ChaUenging for the data analysis are filter channels con- 
taining more than one nitrogen band, as is illustrated in 
Fig. |4l The only channel with a nearly pure contribution of 
a single nitrogen transition is channel 4 which is measuring 
the prominent 2P(0,0) band at 337.1 nm. Therefore, the 
main task for the data analysis is the separation of the dif- 
ferent contributions in the individual filter channels. This 
can be achieved to some extend by a global analysis of the 
complete set of measurements, making use of the physical 
relations between the different nitrogen bands, as discussed 
in Section [2l Technically this was done by a constrained 
X^-minimization of all datasets, where a dataset consists of 
the synchronous data of all filter channels of a run with a 
certain gas at a certain pressure. The most universal con- 
straints directly follow from the general decay law ^ and 
require identical lifetimes Ty/ (p, T) and intrinsic yields Y^, 
for nitrogen bands emerging from the same initial state v' 
in all filter channels of a dataset. Furthermore, the inten- 
sity ratios Rv',v" of transitions from the same initial state 
v' must have the same values in all datasets and therefore 



for each nitrogen band in the different filter channels, have 
been calculated assuming a line spectrum and are listed in 
Table [1] The corresponding systematic errors are the ma- 
jor source of uncertainty for the absolute determination of 
the fluorescence yield. Further details about the experiment 
and the calibration procedures can be found in [9] . 

4. Measurements &: Data Analysis 

The data presented in this work were taken from August 
to November 2005 consisting of about 50 runs in pure nitro- 
gen (/n2 = 1.00), dry artificial air (/nj = 0.78, /02 = 0.21, 
/ai- = 0.01) and a nitrogen-oxygen mixture (/n2 = 0.90, 
/02 = 0.10) to reveal the quenching effects of the different 
air constituents. In addition several runs with pure nitro- 
gen plus a variable amount of water vapor have been per- 
formed. The measurements were done at pressures ranging 
from 2 hPa to 990 hPa at room temperatures between 15 °C 
and 23 °C, varying with the outside weather conditions. In 
order to accumulate sufficient statistics, a single run lasted 
between 12 and 30 hours, depending on the gas mixture and 
the pressure. Before each run the chamber was flushed and 
filled with fresh gas to avoid aging and contamination ef- 
fects. Extensive calibration measurements jSj were carried 
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link different runs with each other. These constraints are 
always valid since they do not depend on any assumption 
about the excitation or de-excitation mechanism. To fur- 
ther improve the results additional, model dependent, con- 
straints may be applied such as the linear pressure depen- 
dence of the reciprocal lifetimes as follows from Eq. lfTO|l . 
This global approach directly results in a consistent de- 
scription of the fluorescence yield with a minimal set of 
parameters. Accordingly, the global minimization function 
was constructed as a sum over all datasets d, filter channels 
c and time bins b of the corresponding time distributions: 



EEE 

d c=l b 



Udcb — {Ode + Bdcb + Sdcb) 



'Vdcb 



(23) 



In this expression, ydcb denotes the bin content with statis- 
tical error Oy^^^^ = -Jyld, and Ode corresponds to the con- 
stant noise offset per bin of filter channel c in the dataset d. 
The noise offset Ode was individually determined for each 
time distribution (see Fig. [5|) by averaging over the bins 
between 2 ns and 12 ns in front of the rising edge of the flu- 
orescence signal. The electron correlated background Bdcb 
was determined and parametrized by means of background 
measurements with a completely evacuated chamber. Its 
origin was not fully understood, however it is only relevant 
for low pressure measurements [9j. Finally the total fluo- 
rescence signal Sdcb in each bin & is a sum over all nitrogen 
transitions contributing to this channel 



f^d.ch — t n 



, \ ^ cv'v 

^dcb — ^DAQ ' 2^ ^det 



^b {Rv' ,v" ■Nd Tdv') 



(24) 



where the summation is done over all initial electronic- 
vibrational states v' as well as over the final states v". Note 
that in the case of contributions from different electronic 
band systems the summation over v' includes the summa- 
tion over the electronic band systems as well! The function 
^b {Rv' ,v" • Ndv' , Tdv' ) is the integral of Eq. (|17p over the 
width of bin b and depends on the lifetime Tdv' of the initial 
state v' and on the total number of emitted photons by the 
transition v' — *■ v" which are related to the total amount 
of photons Ndv' of the main-transition by Rv' ^v" • Ndv'- In 
addition, even if not explicitly quoted in Eq. (IMl) , F\j de- 
pends on the time origins tg and time resolutions at of each 
channel which are kept at the same values in all datasets 
during the minimization procedure. In order to retrieve the 
number of sampled signals Sdcb, each integral function Fb 
must be multiplied by the detection efficiency e'^^^^ of the 
corresponding transition and the data acquisition efficiency 
^'daq (duty cycle) of the filter channel. The detection ef- 
ficiencies have been determined by calibration measure- 
ments, as explained in the previous section, and are listed 
in Table [1] whereas the data acquisition efficiencies s^jaq 
are simultaneously measured during a run using the scaler 
values and were on the order of 91 % for pure nitrogen and 
96 % for dry air 0. 



The only parameters which are finally varied during the 
global minimization procedure thus are the total number of 
photons of the main-transitions Ndv' , the intensity ratios 
Rv',v" and the lifetimes Tdv> of the initial states v' . If the lin- 
ear pressure dependence of the reciprocal lifetimes is taken 
as an additional constraint the lifetime Tdv' in Eq. (|24j) is 
substituted by relation ^ and the intrinsic lifetime Tq as 
well as the various quenching rate constants Q]^ (T) are var- 
ied instead. In this way the consistency of the parameters 
is always ensured. 

According to the separation of the excitation and de- 
excitation processes in Eq. the following data analysis 
was subdivided into two steps. First, the de-excitation 
parameters such as intrinsic lifetimes Tq , intensity ratios 
Rv' .v" as well as the various quenching rate constants 
(T) for the different air constituents, were determined 
using the full usable energy range between 250 keV and 
2000 keV to maximize statistics. Afterwards, the obtained 
values for the de-excitation parameters were fixed and 
the intrinsic fiuorescence yields Y^, {E) of the correspond- 
ing main-transitions were determined in different energy 
intervals, to study their energy dependence. 



4.2. Study oj De-Excitation Parameters 

The present study concentrates on the eight strongest 
nitrogen transitions listed in Table [1] which belong to 
the three independent vibrational sub-systems 2P(0,w"), 
2P(l,t;") and lN(G,w"). To obtain stable and rehable results 
by the global minimization procedure explained above, at 
least two nitrogen transitions from different vibrational 
sub-systems need to be measured separately. With the 
current filter set this was only possible for the 2P(0,0) 
transition in channel 4. Nevertheless, the 2P(1,0) transi- 
tion in channel 1 can also be approximately considered as 
a pure transition since it contains only small contributions 
of the 2P(2,1) and 2P(3,2) bands, as can be seen in Fig.H 

An additional challenge was the measurement of the 
1N(0,0) transition in channel 5, since it is contaminated 
by several weak bands of the higher vibrational transitions 
2P(3,6), 2P(2,5) and 2P(1,4). Fortunately, transitions of 
the IN system have much longer intrinsic lifetimes and are 
also more affected by the coUisional quenching than the 2P 
transitions. Therefore the evolution of their time spectra 
with pressure differs much from the 2P transitions, what 
helps to separate both band systems from each other. Since 
only 2P transitions with u' = 0, 1 are considered in this 
analysis, all the other 2P contaminations in channel 5 with 
v' > 1 are merged into the 2P(1,4) transition assuming 
their pressure dependence to be similar. In this way, the 
2P(1,4) transition becomes an effective measure for all 2P 
transitions in channel 5 and is expected to be more intense 
than the pure 2P(1,4) transition. In order to avoid negative 
effects on the results for the 2P transitions in the other 
filter channels, the analysis was done in two steps. At first, 
only the 2P-transitions in channels 1 to 4 were simulta- 
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(a) Nitrogen at 800 hPa (smaller time scale!) 
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Fig. 7. Time spectra in pure nitrogen at 20 hPa and 800 hPa (different time ranges are displayed in both cases). The signal in channel 
is not fitted but has been reconstructed from the signals in the other channels (see also Fig. If OH . For low pressures, the electron correlated 
background contributions (pink line) are not negligible and have been parametrized by an empirical formula 9 . 




Fig. 8. Stern- VoUmer plots of the reciprocal lifetimes for the vibrational states v' = 0,1 of the 2P system. The lines are not a direct fit to 
the data points but illustrate the pressure dependence derived from the lifetimes and the quenching rate constants, which were obtained by 
a global minimization with a linear constraint on the reciprocal lifetimes. 



neously studied by excluding channel 5 from the analysis. 
Afterwards, just channel 5 was analyzed, keeping the life- 
time of the 2P(1,4) transition in each dataset fixed to the 
values obtained by the previous analysis of the 2P system. 

An example for the resulting contributions in each filter 
channel is given in Fig. [7] where the time spectra of two 
datasets in pure nitrogen at pressures of 20 hPa and 800 hPa 
are shown. Note that for high pressures the time scale is 
much smaller than for the low pressure measurements. At 
20 hPa the 1N(0,0) band in channel 5 is clearly dominat- 
ing whereas at high pressure the 2P contaminations are of 
the same order than the 1N(0,0) signal. Channel with 
the M-UG6 broad band filter is always excluded from the 
minimization procedure since it may contain other contri- 
butions from unknown or disregarded fluorescence transi- 
tions. Nevertheless, its signal can be reconstructed with 
the information from the other channels providing a useful 
crosscheck of the minimization results as is shown in Sec- 
tion 

4.2.1. Analysis of the 2P system 

As explained above, only the filter channels 1 to 4 were 
used to study the quenching parameters of the 2P system. 
Furthermore, the analysis was restricted to measurements 
at pressures higher or equal than 10 hPa to avoid system- 
atic effects due to slight gas impurities. The time spec- 
tra of 46 runs in nitrogen, artificial air, and the nitrogen- 
oxygen mixture have been simultaneously fitted by the gen- 
eral minimization function (|23p . In the first instance the 
minimization function was only constrained by the univer- 
sal relations following from the decay law (jH). The result- 
ing lifetimes for the vibrational levels v' = 0, 1 are recip- 
rocally plotted versus the pressure in the Stern- VoUmer 
plots shown in Fig. [5] For all gases the data points show a 
clear linear dependence on the pressure as expected from 
Eq. (fT0|) . Some outliers may be explained by an insufficient 



separation of different contributions by the minimization 
procedure. Such effects can be reduced by further constrain- 
ing the lifetime in Eq. ([24|) through the linear expression ([9]), 
which additionally allows to account for slight temperature 
variations between different runs. The results of a second 
minimization, using these constraints, are represented by 
the lines drawn in Fig.[8l The y-axis intercept is in common 
for all lines within a graph, indicating the reciprocal intrin- 
sic lifetime of the corresponding initial state v'. The lines 
are in good agreement with the single data points obtained 
by the previous minimization without the linear constraint. 

This final minimization, with a global reduced Xr of 1-77, 
leads to a consistent set of intrinsic lifetimes Tq , quench- 
ing rate constants and intensity ratios R^' ,v" which are 
listed in Tables [2] - 14] together with the results of other au- 
thors. The value of 32.9 ns for the intrinsic lifetime of the 
2P(l,u") sub-system seems to be too low compared to the 
other authors. The reason might be the contamination of 
the 2P(1,0) main-transition in channel 1 with slight contri- 
butions of the 2P(2,1) and 2P(3,2) transitions which can- 
not be separated by this analysis. This contamination is 
also supposed to be the main reason for the deviation of 
the intensity ratios of the 2P(l,ti") sub-system from the 
theoretical expectations of Gilmore et al. in Table U) 

4.2.2. Analysis of the IN system 

The 1N(0,0) transition is the only strong transition of the 
IN system between 300 nm and 400 nm which was stud- 
ied in this work. Due to the reasons explained before, the 
quenching parameters of this transition were determined 
by a separate analysis of channel 5 only. To separate the 
1N(0,0) transition from the 2P contaminations, the life- 
time of the 2P(1,4) transition was fixed in each dataset to 
the corresponding value taken from the previous analysis 
of the 2P system, and only the intensity ratio i?i.4 was al- 
lowed to vary during the global minimization procedure. 
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Fig. 9. Stern- VoUmer plot of the reciprocal lifetimes for the lN(0,v") 
sub-system. As in Fig. |8] the lines correspond to the pressure de- 
pendence obtained by a global minimization with a linear constraint 
on the reciprocal lifetimes. The deviation between the single data 
points and the minimization result is due to unseparated 2P-con- 
taminations. 

In contrast to the 2P system, the pressure range could be 
extended down to 2 hPa since the IN transitions become 
much stronger at lower pressures and are also less affected 
by gas impurities due to the large self-quenching rate of ni- 
trogen. As for the 2P system the minimization was done 
in two steps with increasing levels of constriction. The first 
minimization, using only the most general constraints, re- 
sults in the individual data points for the intrinsic lifetimes 
shown in the Stern- VoUmer plot in Fig. [9] Compared to 
the corresponding plots of the 2P system in Fig. [3 a linear 
behavior is hardly visible. This is caused by an insufhcient 
separation of the 2P contaminations, for pressures larger 
than 50 hPa. Below this point the data seem to adopt a 
linear behavior since the intensity of the 1N(0,0) transi- 
tion becomes strong enough to be silhouetted against the 
2P contaminations, as can also be seen in the time distri- 
butions of Fig. [71 Therefore only runs between 2 hPa and 
50 hPa were considered in the second minimization step 
in which the lifetimes were additionally constrained by the 
linear expression as was already done for the 2P sys- 
tem. This minimization resulted in a global reduced Xr of 
1.26 and a set of de-excitation parameters which are listed 
in Tables [5]- 13] together with the values for the 2P system. 
It must be stressed again that the intensity ratio of the 
2P(1,4) transition in Table |4] is just an effective measure for 
all 2P contaminations in channel 5. Therefore its value is 
expected to be larger than for the pure 2P(1,4) transition. 

The results of this final minimization are represented by 
the lines in Fig. [51 which are clearly deviating from the sin- 
gle data points obtained by the first minimization without 
the linear constraint. Therefore, the separation of the 2P 
from the IN contributions in channel 5 only seems to be 
possible by using all the known physical relations to con- 
strain the minimization function. Furthermore, nearly no 
difference between the three gas mixtures can be recog- 
nized, since the slopes of the lines turn out to be roughly 



Fig. 10. Relative differences between the reconstructed and the mea- 
sured signal in channel with the M-UG6 filter (systematic errors 
are included). From this follows that the contribution of neglected 
nitrogen bands is less than 4 %. 

the same. This leads to the conclusion that the quenching 
of the lN(0,u") sub-system is mainly due to nitrogen self- 
qucnching, since the different contributions of oxygen show 
hardly any effect on the lifetimes. 

4.2.3. Comparison to M-UG6 Data 

The de-excitation parameters determined in the last two 
sections completely describe the pressure dependence of 
the relative intensities for the investigated nitrogen transi- 
tions. To verify their applicability on the whole atmospheric 
pressure range, the time spectra of channels 1 to 5 of the 
complete set of measurements were fitted again with the 
de-excitation parameters fixed to the values listed in Ta- 
bles [5]- [31 Only the intensities of the main-transitions were 
allowed to vary freely in each dataset. Afterwards, the flu- 
orescence signal in channel with the M-UG6 broad band 
filter was calculated using the fitted intensities and the de- 
tection efficiencies (see Table [Ij of each transition. An ex- 
ample for the reconstructed time distributions in channel 
is given by the upper left graphs of Fig. [71 for pure nitrogen 
at two different pressures. The yellow histograms represent 
the measured time distributions and the calculated contri- 
butions of each nitrogen transition, as well as their sum, are 
plotted on top of the measured distribution. For both pres- 
sures, the signal in channel could be almost completely 
recovered by the information from the other channels. This 
worked also for the other gas mixtures and pressures, as is 
shown in Fig. llGl where the relative differences between the 
reconstructed and measured integrated fluorescence signals 
in channel are plotted versus the pressure. The systematic 
uncertainties are already included in the error bars and the 
large fluctuations of the data points below 10 hPa can be 
explained by the limited statistics for low pressure runs. For 
pressures larger than 10 hPa the data points stabilize and 
the reconstructed signals are about 4 % smaller, but still 
compatible with the measured signals. However, a small ex- 
cess of the measured signal in channel is expected, since 
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Fig. 11. Reciprocal lifetimes in humid nitrogen versus the water vapor partial pressure. The nitrogen partial pressure was always kept at 
30 hPa. The data points at PH2O = correspond to the reciprocal lifetimes l/ry in dry nitrogen at a pressure of 30 hPa. 



this channel is also measuring transitions which were ne- 
glected in the present analysis. From this follows that the 
contribution of disregarded transitions to the total nitro- 
gen fluorescence spectrum in the M-UG6 filter range is less 
than 4 %. 

4.2.4. The Effect of Water Vapor 

The influence of water vapor on the nitrogen quench- 
ing was studied by adding different concentrations of 
water vapor to pure nitrogen at a (partial-) pressure of 
PN2 = 30 hPa. The water vapor concentration was mon- 
itored by a calibrated humidity probeff] which was mea- 
suring the relative humidity with an absolute systematic 
uncertainty of 1.5 % rH. The water vapor partial pressure 
PH2O was derived from the relative humidity according 
to the magnus formula 19116] . Six runs were performed at 
room temperatures between 15 °C and 17 °C and relative 
humidities ranging from 10 % rH to 60 % rH. Each run 
lasted for roughly 30 hours but only data of the last ~ 20 
hours was used for the analysis, in order to ensure the wa- 
ter vapor concentration to be in thermal equilibrium with 
the surrounding chamber walls. 

As before, the data analysis was done by simultaneously 
fitting the time spectra of all datasets with the minimiza- 
tion function (|23p , by applying only the most general con- 
straints and by fixing the intensity ratios to the values ob- 
tained from the previous analysis listed in Table U) The 
global reduced Xr of this minimization was 1.09 and the 
resulting reciprocal lifetimes for the three vibrational sub- 
systems are drawn versus the water vapor partial pressure 
PH2O in Fig- [HI The data points at PH2O — correspond 
to the results of the previous analysis for dry nitrogen at a 
pressure of 30 hPa. According to Eq. ^ the water vapor 
quenching rate constants Qh^q were determined through a 
fit of the linear expression 



PN2 
kT 



PH2O 
" kT 



Qh2 



(25) 



= : 1/fo 



to the data points in Fig. [TI] resulting in the values listed 
in Table [3] For the 2P system the water vapor quench- 
ing rate constants were roughly two times larger than the 
corresponding values for oxygen, whereas the quenching 
of the IN system by water vapor is roughly three times 
stronger than by nitrogen or oxygen. This large value may 
be explained by the polar nature of the H2O and the 
molecules. Despite of the large values for the water vapor 
quenching rate constants the effect on the nitrogen quench- 
ing in the atmosphere is expected to be rather small since 
the water vapor concentration in air usually is much lower 
than the nitrogen or oxygen concentrations. 

4.3. Study of Energy Dependence 

According to Eq. (fT6|l the total number of photons N^^k, 
emitted by a main-transitions v' — s- fc, is related to the 
intrinsic fluorescence yield Y°, ^ {E) through the expression 



N„,,k{p,T,E) 



Tv' {P, T) 



Y:,[,{E) ■ E,,,{E) 



(26) 
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The lifetime t^' {p, T) is only affected by the de-excitation 
parameters, which are known from the analysis described 
before. Therefore the last missing quantity for the deter- 
mination of the intrinsic fluorescence yield is the deposited 
energy Edep{E) in the chamber volume which was obtained 
from Monte-Carlo studies described below. 

4.3.1. Energy deposit in Chamber 

A precise knowledge of the energy deposit in the cham- 
ber is crucial for the correct determination of the fluores- 
cence yield. In confined volumes, the energy deposit is not 
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Table 2 

Intrinsic lifetimes Tq in [ns] for the initial states of the investigated vibrational sub-systems. Only 
statistical errors are quoted. 



2P{v' = 0) 


2P(v' = 1) 


1N{»;' = 0) 


Reference 


38.9 ±0.3 


32.9 ± 0.5 


65.2 ± 18.7 


this work 


42.0 ±2.0 


41.0 ±3.0 




Pancheshnyi et al. [17) 


41.7 ± 1.4 


41.7±2.1 




Morozov et al. [18] 


44.5 ±6.0 


49.3 


65.8 ±3.5 


Bunner [3] 


37.1* 


37.5* 


62.3* 


Gilmore et al. |14| 



* Theoretical values 



Table 3 



Quenching rate constants QJ^. in [10 
sub-systems, given for T = 293K {~ 


^"cm^s ^] for the initial states of the investigated vibrational 
20° C). Only statistical errors are quoted. 


Molecule 


2P{v' = 0) 


2P{v' = 1) 


1N{?;' = 0) 


Reference 




0.11 ± 0.00 


0.29 ± 0.00 


5.00 ±0.17 


this work 




0.13 ±0.02 


0.29 ±0.03 




Pancheshnyi et al. |17| 




0.12 ± 0.01 


0.25 ± 0.01 




Morozov et al. [TS] 




0.10 ±0.01 


0.22 ± 0.03 


3.80 ±0.26 


Nagano et al. [5] 




0.07 


0.23 


2.90 


Bunner 


O2 


2.76 ±0.01 


2.70 ±0.03 


5.24 ±0.79 


this work 




3.00 ±0.30 


3.10 ±0.30 




Pancheshnyi et al. 1171 




2.62 ± 0.19 


2.77 ±0.45 


2.39 ±0.40 


Nagano et al. [5] 




1.35 




8.38 


Bunner '3' 


H2O 


5.43 ± 0.12 


5.78 ± 0.17 


16.02 ± 1.09 


this work 




3.90 ±0.40 


3.70 ±0.40 




Pancheshnyi et al. |17j 




7.10 ±0.70 


6.70 ±0.70 




Morozov et al. 18 



Table 4 

Intensity ratios ^// of the investigated transitions. The values of Bunner have been derived from the 
fluorescence efficiencies quoted in his work. The second error corresponds to systematic uncertainties 
of the detection efficiencies. 



Transition 


this work 


Bunner [3I 


Pons et al. [19] 


Gilmore et al. [Ti] 


2P(0,0)'' 


1.00 ±0.00 ±0.00 


1.00 


1.00 


1.00* 


2P(0,1) 


0.69 ±0.00 ±0.06 


0.75 


0.63 


0.67* 


2P(0,2) 


0.29 ± 0.00 ± 0.02 


0.26 


0.25 


0.27* 


2P(1,0)" 


1.00 ±0.00 ±0.00 


1.00 


1.00 


1.00* 


2P(1,2) 


(0.33 ± 0.03 ± 0.04)'' 


0.58 


0.45 


0.47* 


2P(1,3) 


(0.34 ±0.01 ±0.03)'' 


0.54 


0.43 


0.41* 


2P(1,4) 


(0.46 ± 0.02 ± 0.04)^= 


0.32 


0.16 


0.20* 



* Theoretical values. 

" Main-transitions (Ratio equals to one per definition). 
'' Biased value due to contamination of 2P(1,0) band. 
Biased value due to contamination of 2P(1,0) and 2P(1,4) bands. 
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Fig. 12. GEANT4 simulation of the spatial distribution of the de- 
posited energy density within the AirLight chamber at a pressure 
of 500 hPa, for detected electron energies between 250 keV and 
2000 keV. 

necessarily equal to the ionization energy loss, as it follows 
from the Bethe-Bloch formula, since secondary delta elec- 
trons may carry away part of the energy. Furthermore, ad- 
ditional secondary electrons emerging from collisions of the 
primary electrons with the collimator walls, or electrons 
repulsing out of the scintillator, may even increase the ef- 
fective energy deposit. Therefore the energy deposit in the 
chamber is expected to deviate from the ionization energy 
loss. Energy losses due to Brcmsstrahhmg in the gas are 
negligible at these energies. 

To determine the mean energy deposit in the chamber 
of an electron which was detected with a certain energy 
in the scintillator, detailed GEANT4 simulations have 
been performed 9]. The simulations were done using 
GEANT4 7.1 20J and the low energy extension for the 
electromagnetic processes. The simulated ^°Sr electron 
energy spectra were convoluted with the energy resolution 
function of the scintillation detector and compared to the 
measured energy spectra. The "facrange" parameter of 
the multiple scattering model, which affects the penetra- 
tion depth of the electrons into a new material, had to be 
tuned to a value of 0.01. This significantly improved the 
treatment of the electron backscattering within the lead 
collimator which has to be taken into account in order to 
produce realistic energy spectra [9 . After the simulation 
was tested and verified with real electron spectra, several 
simulations for dry air, at all pressures applied to the real 
measurements, were performed. 

The simulation directly generated two-dimensional his- 
tograms of the deposited energy in each volume element of 
the chamber, with a spatial resolution of 1 mm^. Each elec- 
tron, no matter if it is the primary or a secondary electron, 
contributes to the energy deposition. An example for such 
an energy deposit map is given in Fig. [T^] for detected elec- 
tron energies ranging from 250 keV to 2000 keV. If the num- 
ber of emitted fluorescence photons is proportional to the 
deposited energy. Fig. [12] can be reinterpreted as the spa- 
tial distribution of the fluorescence emissions in the cham- 



Fig. 13. Mean deposited energy per electron in the chamber versus 
the detected energy in the scintillator for various pressures. The red 
lines correspond to the ionization energy loss in dry air given by the 
Seltzer &; Berger formula [21) . 

ber. This was done to determine the angular acceptance £n 
of the PMTs 9 which was used to calculate the detection 
efficiencies in Table [T] 

In order to study the energy dependence of the energy de- 
position, several energy deposit maps for successive, 25 keV 
wide energy intervals were generated. Each map was inte- 
grated over the volume of the chamber and divided by the 
number of detected electrons to derive the mean energy de- 
posit in the chamber per detected electron. In Fig. fT51the 
mean energy deposit is plotted for several pressures ver- 
sus the detected electron energy in the scintillator. The 
comparison to the Seltzer & Berger description 2l\ for the 
mean ionization energy loss reveals a 15 % to 20 % decrease 
for electrons detected with energies larger than 500 keV. 
Below this point the deposited energy seems to agree with 
the ionization energy loss suggesting most of the secondary 
electrons to be stopped within the gas volume of the cham- 
ber. 

4.3.2. Intrinsic Yields 

To determine the intrinsic yields Y^, ,^{E) and to study 
their energy dependence the measured data were divided 
into seven sub-samples of 250 keV energy intervals covering 
the range between 250 keV and 2000 keV. Each sub-sample 
was fitted separately according to the minimization proce- 
dure explained in Section [4. II During the minimization all 
the de-excitation parameters were fixed to the values listed 
in Tables [2] -|4] and only the total number of photons Njj>^k 
emitted by the main-transitions was allowed to vary freely. 
Afterwards the intrinsic fluorescence yields Y^, ^ (E) in each 
energy interval were evaluated by solving Eq. (|26p using 
the total number of photons Nyi^k from the fit and the de- 
posited energy, which was calculated from the total num- 
ber of detected electrons multiplied by the mean deposited 
energy per electron from the simulation studies explained 
above. 

An example for a medium pressure of 400 hPa is given in 
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Fig. 14. Energy dependence of the intrinsic fluorescence yields at 
a pressure of 400 hPa. The horizontal lines are not a direct fit to 
the data points, but represent the minimization results using the 
complete energy range. 



Fig. 15. Pressure dependence of the intrinsic fluorescence yield which 
presumably is due to experimental deficiencies. The horizontal lines 
correspond to the values quoted in Table [6] and denote the weighted 
average of the data points. 



Fig. where the intrinsic fluorescence yields of the three 
main-transitions arc plotted versus the electron energies. 
The relative differences between the intrinsic yields in dif- 
ferent nitrogen mixtures seem roughly to scale with the 
different nitrogen concentrations. Furthermore the intrin- 
sic fluorescence yield appears to be constant over the whole 
energy range which has been studied. Some fluctuations at 
high energies can be explained by the reduced statistics due 
to the /3-shape of the electron energy spectrum. In partic- 
ular between 250 keV and 1000 keV, where the ionization 
energy loss varies about 30 %, there is no notable change of 
the intrinsic fluorescence yields. Therefore, the above anal- 
ysis was repeated using the whole usable energy range be- 
tween 250 keV and 2000 keV. The results are represented 
by the horizontal lines in Fig. [14] and agree well with the 
corresponding data points of the single energy intervals. 

To verify the results, this procedure was repeated for 
several pressures even if a pressure dependence was not ex- 
pected from the model described in Section[2l It turned out 
that the intrinsic fluorescence yield apparently decreases 
for lower pressures as is shown in Fig. 1151 The effect is 
strongest for pure nitrogen and seems to decrease with in- 
creasing oxygen concentrations in the other gas mixtures. 
It could not be clarified if this behavior is a real physical 
effect or if it is due to experimental issues. Therefore, the 
final intrinsic fluorescence yield values in Table[6]have been 
calculated as the weighted average over the data points at 
different pressures, and are represented by the horizontal 
lines in Fig. 1151 The corresponding systematic uncertain- 



Table 5 

Systematic uncertainties due to the observed pressure dependence 
of the intrinsic fluorescence yield. 



Gas 


2P(0,0) 


2P(1,0) 


1N(0,0) 


N2 


12 % 


4 % 


5 % 


N2 : O2 (90:10) 


5 % 


3 % 


5 % 


Dry Air 


7 % 


2 % 


4 % 



ties have been determined as the standard deviation of the 
single data points to the averaged value and are listed in 
TableO These uncertainties have been quadratically added 
to the uncertainties of the detection efficiencies in Table [T] 
to determine the systematic errors quoted in Table El Ac- 
cordingly the final systematic uncertainties for the individ- 
ual nitrogen bands in air settle between 14 % and 16 %. 

5. Results and Discussion 

At the beginning of this paper an alternative approach to 
parametrize the nitrogen fluorescence process in air was in- 
troduced. This approach accounts for all physical relations 
between the different nitrogen emission bands and consis- 
tently describes the whole process with a minimal set of 
parameters. The excitation and de-excitation processes are 
completely separated by this model and the correspond- 
ing parameters have a clear physical meaning. The fluores- 
cence yield in this work is defined in terms of photons per 
deposited energy. This definition is much closer related to 
the experimental quantities than the traditional definition 
in terms of photons per meter track length and addition- 
ally allows to account for effects of escaping delta electrons 
in small volumes. The latter issue especially is important 
for the absolute determination of the fluorescence yield in 
small laboratory experiments. 

The AirLight experiment synchronously measured the 
eight strongest nitrogen emission bands between 300 nm 
and 400 nm for pressures ranging from 2 hPa to 990 hPa 
with a systematic uncertainty of about 15 %. Furthermore 
the influence of water vapor on the nitrogen quenching has 
been studied. The experimental data was analyzed accord- 
ing to the model mentioned before. In the investigated en- 
ergy range between 250 keV to 2000 keV the fluorescence 
yield was shown to be independent from the energy of the 
ionizing electrons. 
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Table 6 



Absolute values of the intrinsic fluorescence yields Y^, ^ in [pliotons/keV] for the investigated main-transitions v' ^ k. The values of this 
work are valid for electron energies ranging from 250 keV to 2000 keV. The values of Nagano et al. have been converted from his ^''-values 
according to = The second error corresponds to systematic uncertainties. 



Gas 


2P(0,0) 


2P(1,0) 


1N(0,0) 


Reference 


Na 


0.513 ±0.001 ±0.096 


0.236 ±0.001 ±0.034 


1.513 ±0.014 ±0.223 


this work 




0.318 ±0.005 ±0.041 


0.129 ±0.007 ±0.017 


0.397 ±0.028 ±0.052 


Nagano et al. [6] 


N2 : O2 (90 : 10) 


0.429 ±0.001 ±0.065 


0.210 ±0.001 ±0.030 


1.360 ± 0.009 ± 0.200 


this work 
Nagano et al. 


Dry air 


0.384 ±0.001 ±0.061 


0.195 ±0.001 ±0.028 


1.190 ±0.008 ±0.171 


this work 




0.272 ± 0.007 ± 0.035 


0.122 ±0.007 ±0.016 


0.303 ±0.012 ±0.039 


Nagano et al. |6] 



5.1. Usage of the Model 



To calculate the fluorescence yield Yv>y>{p, T) for a cer- 



tain nitrogen transition v' - 
and pressure p the formula 



at a given temperature T 



(27) 



applies, where the value of the associated intrinsic fluores- 
cence yield Y^, is listed in Table [5] and the values for the 
intensity ratio R^' y and the intrinsic lifetime ro„' are tab- 
ulated in Table [4] and Table [2] compared to the values of 
other authors. In this context it should also be referred to 
the results of the AIRFLY experiment [7j which have re- 
cently been published. However, it is strongly advised not 
to mix the intensity ratios and intrinsic yields from different 
authors since both quantities are correlated to each other. 
Even if some intensity ratios of this work, especially for the 
2P(l,z;") sub-system, deviate from the values of other au- 
thors, they still yield reasonable absolute intensities if they 
are used together with their associated intrinsic yield val- 
ues. 

The pressure and temperature dependence of the life- 
time Tt,/ (p, T) in Eq. (PT)) can be parametrized according to 
Eq. ^ and result in the following expression for the recip- 
rocal lifetime in air: 



1 _ 1 



P 

kT 



(0.78 •gj^^(T) + 0.21. g^jT)) . 



1 - 



P 



P 



(28) 



This expression already accounts for the water vapor 
quenching if the water vapor partial pressure Pjj^q is 
given, which derives from the relative humidity by means 
of the magnus formula [9ll6j . The temperature dependent 
quenching rate constants Q'^ (T) can be calculated from 
the tabulated values in Table [3] by means of Eq. ([7]). 

On the basis of Eq. ^7} fluorescence yield spectra can be 
calculated for the whole pressure and temperature range 
applying to cosmic ray air shower experiments. An exam- 



ple is given in Fig. [16] where the spectrum for dry air at 
4 km height above sea level is compared to the correspond- 
ing spectrum of Nagano et al. j6j . The pressure and tem- 
perature at 4 km height are computed according to the US 
Standard Atmosphere [55^ and roughly describe the atmo- 
spheric conditions as they appear at the shower maximum 
where most of the energy is deposited in the air. 

At a flrst glance both spectra seem to agree within their 
errors although there are big differences between the intrin- 
sic yield values in Table [51 This is due to the reason that 
the de-excitation parameters (intrinsic lifetimes, quench- 
ing rate constants) of both experiment are not exactly the 
same and the intrinsic yield values are calculated in the 
limit of zero pressure, in order to separate the excitation 
and de-excitation processes. In pressure ranges where the 
actual measurements have been performed the agreement 
between both experiments is much better since the differ- 
ent de-excitation parameters compensate for the differing 
intrinsic yield values. 

The biggest difference in Fig. [THj occurs for the 2P(1,4) 
transition, which in the present work appears to be too 
large. This was expected since this work's 2P(1,4) transi- 
tion is still contaminated by neighboring bands and there- 
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Fig. 16. Absolute fluorescence yield spectrum of the investigated 
nitrogen transitions in dry air, at atmospheric conditions as they 
appear at shower maximum, compared to the values of Nagano et 
al. 6 . The 2P(1,4) and the 2P(1,0) transitions of this work still 
contain contributions of neighboring bands (see text). 
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Fig. 17. Fluorescence yield of the three main-transitions and the sum 
of all investigated transitions versus height a.s.l.. The lower panel 
shows the relative differences between the values of Nagano et al. 
and this work in more detail. 

fore is only an effective measure of all the 2P contributions 
in the corresponding filter channel. Also the 2P(f ,0) tran- 
sition is known to be slightly contaminated by neighboring 
bands, however, there seems to be no remarkable difference 
compared to the spectrum of Nagano et al.. 

The differences between both spectra become more ev- 
ident if the fluorescence yield values are compared over a 
large atmospheric range. In FiglT7]the fluorescence yield of 
the 2P(0,0), 2P(0,I) and 1N(0,0) main-transitions as well 
as the sum of all 8 investigated transitions is drawn versus 
the height above sea level. In the lower regions of the at- 
mosphere this work's fluorescence yield values for the 2P 
transitions in average are about 10 % larger, whereas the 
value of the 1N(0,0) transition appears to be 6 % smaller 
than the corresponding values of Nagano et al.. These de- 
viations are still within the systematic accuracies of both 
experiments. Above 10 km height the influence of the differ- 
ent de-excitation parameters becomes visible, leading to a 
slower increase of the values of Nagano et al. . Since the flu- 
orescence light emission in EAS occurs mainly in the lower 
atmosphere a 10 % difference between both fluorescence 
models can be expected. 

Water vapor in the atmosphere leads to a reduction of the 
fluorescence yield as shown in Fig. 1181 where this work's flu- 
orescence yield is compared for dry and humid atmospheres. 
Realistic atmospheric profiles, obtained from radiosonde 
measurements |23I24| at the site of the Pierre Auger Obser- 
vatory [T] in Malargiie (Argentina), have been used to com- 
pute the fluorescence yield with and without taking into 
account the water vapor concentration. The solid lines cor- 
respond to the average of 187 radiosonde launches accom- 
plished in all seasons. The dashed lines indicate the upper 
and lower bounds of one standard deviation to the aver- 
age due to the varying water vapor content in the air. At 
ground level the fluorescence yield of the 2P transitions is 
2 % to 5 % lower if the water vapor quenching is taken into 
account. The IN transition is only reduced by 1 % to 2 % 
since its quenching rate in dry air is already relatively large. 



Fig. 18. Relative differences between this work's fluorescence yield 
with and without taking into account the water vapor content 
of the atmosphere. Realistic atmospheric profiles from radiosonde 
launches I23I24I have been used for this comparison. 

These differences are still concealed by the systematic 
uncertainties of the fluorescence yield but might become 
an issue if the fluorescence yield is known with a better 
precision. Currently the main source of uncertainty are the 
quantum- and collection efhciencies of the PMTs which 
have only been estimated according to the manufacturers 
specifications. Therefore, the ongoing work is concentrat- 
ing on the improvement of the absolute calibration of the 
experiment using Rayleigh-scattering of a pulsed 337 nm 
laser, which substitutes the electron beam in the chamber 
^25j . With this end-to-end calibration of the whole setup, 
the absolute uncertainties of the fluorescence yield are ex- 
pected to drop below 10 %. 
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